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Abstract 

Recently Eling and Oz (EO) proposed a formula for the holographic bulk viscosity, 
in ^rXiv: 1103. 16571 , derived from the null horizon focusing equation. This formula 



seems different from that obtained earlier by Gubser, Pufu and Rocha (GPR) in 



arXiv:0806.04071 calculated from the IR limit of the two-point function of the trace 
of the stress tensor. The two were shown to agree only for some simple scaling cases. 
We point out that the two formulae agree in two non-trivial holographic theories de- 
scribing RG flows. The first is the strongly coupled J\f = 2* gauge theory plasma. The 
second is the semi-phenomenological model of Improved Holographic QCD. 
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1 Introduction and Summary 

The bulk viscosity of strongly coupled thermal systems is a quantity of phenomenolog- 
ical importance. On the other hand it is quite diflicult to compute. The main difliculty 
arises from the fact that the bulk viscosity, as one of the plasma deep-infrared trans- 
port coeflicients, is sensitive to the microscopic (ultraviolet) parameters of the theory]^. 
Necessarily, a computation of the bulk viscosity in a given system requires the under- 
standing of its physics over a wide range of scales. It is perhaps not surprising that 

"'^The only exception is a conformal theory, where the scale invariance imposed the bulk viscosity 
to vanish. 
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the first computation of tlie bulk viscosity in gauge tlieory plasmas [|T| was performed 
in the framework of gauge theory/string theory correspondence [0,0]. 

In m it was observed that for a large class of holographic models, the bulk viscosity 
of the strongly coupled plasma satisfied the following bound 



C /I 
> 2 

■ 3 



7] 

where rj is the universal shear viscosity of strongly coupled holographic plasma 
and Cs is the speed of sound waves in plasma. The computation of the bulk viscosity 
which led to ( |1 . 1| ) was based on analyzing the dispersion relation of the sound waves 
in plasma. Alternatively, the bulk viscosity can be computed using the Kubo formula 



4 1 

where Gr is the retarded correlation function of the stress-energy tensor 



:i.2i 



1.3) 



The holographic computations of the correlator (|1.3| ) for a certain class of dual gravi- 
tational models by Gubser, Pufu and Rocha (GPR) was reported in p. It was claimed 
that some of the Einstein-scalar models considered led to a violation of the bound ( |1.1| ). 
On the other hand, in the Improved Holographic QCD model, [1^, the bound (IT 
comfortably obeyed, |^ 



IS 



Recently an alternative expression for the bulk viscosity in strongly coupled plasmas 
with a holographic dual was obtained by Eling and Oz (EO) in ||10|. They have analyzed 
directly the hydrodynamic limit of the equations of motion of a generic Einstein-scalar 
theory and derived a formula for the bulk viscosity that is apparently different from 
the GPR formula. The EO formula is very general and reads 



1-1 V 



ds 



;i-4) 



where i labels different bulk scalars, (f)\ is the value of the i-th scalar at the horizon, and 
are different conserved charged densities. The case comparable with GPR, involves 
a single scalar field and no charge density. In [|1^] the two formulae were shown to give 
the same result in cases where the adiabatic approximation to the equations is valid, 
but their equivalence in more general cases was put in doubt. 
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In this paper we re-analyze the bulk viscosity in two non-trivial holographic theories, 
the bosonic M = 2* theory |rT|| as well as Improved Holographic QCD, |T^,|TB[|. The 



bosonic M = 2* theory is = 4 superYM, with a non-trivial (and equal) mass for 4 of 
the 6 scalars. Improved holographic QCD on the other hand is a semi-phenomenological 
Einstein-scalar theory tuned to match non-supersymmetric Yang Mills theory in the 
large Nc limitQ . 

We point out that the analysis done in 0], when applied to A/" = 2* gauge theory 



plasma [11,14-lS] at high temperatures, agrees with earlier computations reported 
in H, 10, 19|. The agreement is also checked numerically at all temperatures. 

Both in the M = 2* theory and Improved Holographic QCD we confirm the agree- 
ment between GPR and EO formulae for the holographic bulk viscosity. 

2 The GPR formula for the holographic bulk viscosity 

Here we mostly follow |0. For details we refer to the original work. 

Consider a gravitational system, dual to some strongly coupled gauge theory plasma, 
described by an Einstein-dilaton system of the form 



(2.1) 



The black brane background geometry dual to a thermal state of the plasma takes the 
form 

ds^ = gj^Jdz^-dx^ = e2^('-) [-h{r)dt^ + f/f^] + e^^^^)^ , = r . (2.2) 

Notice that the field was chosen as a radial coordinate^. 

^Similar Einstein-scalar theories were also proposed to describe the crossover behavior of QCD 



with light quarks in |20 . 

•^One might worry whether (p is monotonic from the boundary to the black brane horizon. In the 

Einstein-dilaton theory, there are solutions where 4>' vanishes along the flow. These where analyzed 
in and shown to be unphysical, violating the Gubser bound [ p2| . Therefore, this is not expected 
to happen in the middle of an RG flow. At theories with an extra gauge field and at finite density 
however, it is possible that 0' = in a physical solution. A class of such examples were studied 
recently in 1241. 
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The background equations of motion take a simple form 

=A" - A'B' + - , 
6 

=h" + (AA' - B')h' , 

(2 3) 

=6A'h' + h{24A'^ - 1) + 2e^^V , 

h' e^^ 
h h 

To compute the correlation function ( |1.3| ), the authors of considered an S0{3)- 
invariant fluctuation of the metric 5g^y{t, cj)) — )■ e~^^^ g'^^} {(i))H in the gauge 50 = 0. 
It was shown that the equation for Hn decouples from the rest of the fluctuation 
equations and is 

One further has to solve (|2.4|) with the following UV (r — i- 0) and IR (r — )■ 
boundary conditions: 

UV: limi7n = l, (2.5) 

1 — >o 

IR: ifn^cri(0/.-r)-^-/^'^^ + Ox(0;,-r)+^-/^-^, as r ^ 0^ . (2.6) 
The bulk viscosity, computed from (|1.2| ), is given by 

- = . ] lim \c7A\ (2.7) 

where one has to use the universality of the shear viscosity The authors of 

used ( p.3|) to obtain 

and arrived at the final formula for the bulk viscosity ratio 



limlqil'. (2.9) 



^We independently reproduced this equation. We also verified the consistency of the gauge choice 

^Note that in the gauge we are working the position of the black hole horizon is identified with 
the value of the scalar at the horizon 0^. 
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2.1 The EO versus GPR formula for the bulk viscosity 



In Eling and Oz, by analyzing the hydro dynamic limit of the scalar-tensor equa- 
tions, produced the following expression for the holographic bulk viscosityf] 

2 



EO 



ds 



1 



(2.10) 



Even though (up to a factor of Cn) ( |2.10| ) and ( p.7p appear to be the same, they are, 
in fact, different: in ( |2.10| ), 



d (lim^^^^ y4(0)) 



^A'( 



EO 



GPR 



lim 

<t>^4>h 



dA{(P) 



XML 

<A"4 3\/'(0) 



lim 



• (2.11) 



To be specific, in A/" = 2* gauge theory plasma at high temperature (see Appendix 
for some details) 



A'i 



TC 



EO 



mt 



m 



2\ 0^ 



O 



GPR 

From ( p.l2| ) and (|2.13| ) it is clear that ( p.7P would produce the correct expression 
for the M = 2* plasma bulk viscosity, provided[] 



2^2 



2\ 0^ 



(2.12) 



(2.13) 



-111 



2 + 0{ji 



mt 



(2.14) 



Af=2*, prediction 

In the next section we explicitly compute Ci^, and find that it agrees with ( |2.14 ). 



3 Bulk viscosity calculation in the J\f = 2* plasma 

In this section we will first address the calculation of bulk viscosity in the J\f = 2* 
theory. 

^ The formula derived in JTo[ applies also to systems at finite charge density. Here we restrict our 

attention to zero charge density systems. The EO formula was further tested in p5[ . 
^We assume the cj ^ limit taken. 
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3.1 The computation of in x-gauge 

We find it convenient to recast the equation for Hu in terms of x coordinate, defined 
as 

x=l-Vh. (3.1) 

Notice that x — 0+ corresponds to the boundary and a; — t- 1_ to the horizon. 

In this gauge the background equations take the form (all the derivatives are with 
respect to x): 

= A"-4{A'f + -^ + h<P'f, (3.2) 
1 — X o 

where 

The equation for Hu is somewhat complicated 

= H';^ + -Hi H'u + V,2 Hn , (3.5) 

where we collected the coefficients Hi in Appendix 0. In order to compute the bulk 
viscosity (p.7|), we need to solve ( p.5|) subject to the following boundary conditions: 

UV : hm Hn = I , (3.6) 
IR: /Jii ^cri(l-s)-'"/2-^ + X (l-a;)+'"/2-^, as a; ^ 1_ . (3.7) 



For generic w, c^^ 7^ c^^^ (see ( |2.6|) ), however, in the hydrodynamic limit 



limcii = lim qi = • (3.8) 

We can test (p.2|) - (|3.5|) with a simple, exactly solvable background, like the expo- 
nential potential case. This is described in Appendix |C[ 

3.2 M = 2* plasma at high temperatures 

The effective action of the gravitational dual to strongly coupled J\f = 2* plasma with 
a bosonic mass deformation is given by ||18[ 



S=j^ I deV^[\R-^{daf-V], (3.9) 



V = -^e-^" - -e^" . (3.10) 
4 2 ^ ^ 



where the potential is Q 



Notice that the canonically normahzed scalar is = v24a;, and therefore 

K<A=4|^- (3-11) 



We will study the theory ( p.9|) in the high-temperature regime. In this case (see 
Appendix ^ 

= p = 1 + Siai + 0{Sl) , 

1 (3.12) 
A = ln53 - - ln(2a; - x^) + 5?^! + 0(5^) . 

In the hydrodynamic limit, i.e., u ^ 0, and to leading order in Si, we find 

n ^ ^ n" a. (a^«U2 - x)(x^ - 2x + 4) - 2ai(l - x) 
+ 0{6,) =H,, + «;^2(i_^)(2_^)2 



2((4x-2x2)a; + ai(x- 1)) 
^ (2 - x)2(l - x)2«;x2 

Notice that there is dependence only on ai, which satisfied the following equation 

Even though we know an analytic solution for ai (see ( [A.2| )), we can not solve for 
Hii analytically. We find it convenient to use numerical techniques to solve both ( |3.13| ) 
and ( p.l4| ). Near the boundary we have 

/ oo 1 \ 

ai = y/x ^ ^ a„,fc In'' X , (3.15) 

\n=0 A:=0 / 

with normalization^ ao,i = 1, and 

an,fc = an,fc^ao,o^ • (3.16) 

For example, for the first few terms we have: 

_11 _1 _55 _5 

cti.o — 2 + ^'^0,0 , ^1,1 — -, ~ 16 ~'~ 32*^°'° ' '^^'^ ~ 32 ' (3-1' j 



®We set the five-dimensional gauged supergravity coupfing to one. This corresponds to setting the 
radius i of the five-dimensional sphere in the undeformed metric to 2. 

^The overall normalization of ai is arbitrary, we choose the leading In a; coefficient to be 1. 
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The asymptotic expansion for i^n is a bit unusual because the perturbing operator has 
scahng dimension 2: 

oo n+1 ^ 

^11 = E E ^".^ ^^-fTZu: ' (3-18) 



n=0 fc=0 

with normahzation /io,o = 1, see (3^). Here, 



(ao,o + 2 + Inx)'' 



hn,k ^71, fc ( 



'0,1 



(3.19) 



For the first few terms we have: 



hi. 



-ho. 



hi,2 — ~-ho,i 



h2,o 



h2,i = ^(1 - ho,i) , 



^2,2 — 77:^0,1 

lb 



hifi = -1 , 
1 

"4' 

Near the horizon, y = 1 — x, we obtain 

a, = a^f^a^y- = ( 1 - - ^-^ - + 



(3.20) 



^2,3 — ^^0,1 



n=0 

00 



4' 



64^ 



256' 



/ 1 

n=0 ^ 



y 



27 



(3.21) 



64" 1024 
Altogether we have four integration constants: 

{ 0.0,0 ! ^0,1 ! '^'^ , h^ } , 

precisely what is needed to solve uniquely the system of two second order ODEs: ( p.l3|) 
and ( p. 14 ). Using numerical techniques developed in |26] we find 



ao^o = -2.079441(5) , a'' = -2.221441(5) 
/io^ = -2.000000(0) , /i^ = 2.000000(0) . 

Of course, {0,0,0,0,'^} are known analytically from ( |A.2|) , 



(3.22) 



TT 



{ao,o, a^} = {-ln8,--^ } 



and are in excellent agreement with ( p.22|) . 
From (|]22D, 



'11 



h'' 



(3.23) 



Af=2* 



to a very good accuracy, confirming the agreement of bulk viscosity for the high- 
temperature M = 2* plasma from ( p. 71 ) with earlier computations [^|10|,|19 . 
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Figure 1: (Color online) Comparison of the GPR prediction for M = 2* plasma bulk 
viscosity with the explicit computations from the quasinormal modes 0]. The dashed 
vertical green line represents the critical point of the theory ^ = 5c = 5.4098(6) 
associated with the second-order phase transition p7|,p8| . 



3.3 Af = 2* plasma at generic temperatures for ^ > 



It is straightforward to extend the analysis of the previous section to generic temper- 

2 

atures M = 2* gauge theory plasma for physical mass deformations, i.e., > The 
background geometry was studied in [jl8[, and the bulk viscosity (from the sound waves 



dispersion relation) was computed in [Q. The results of the analysis are reported in 
Figure |l|. We further verified that the GPR formula (|2.9| ) for the bulk viscosity, when 
applied to M = 2* plasma, agrees with the bulk viscosity of the theory at criticality @] 
computed from the sound waves dispersion relation to ~ 5 x 10"''. 



4 Bulk viscosity calculation in Improved Holographic QCD 



In this section we perform an independent calculation of the coefficient by the 
methods developed in |^T],^. As described in section 7 of one can work out the 
thermodynamics of gravity- scalar system entirely by solving a system of coupled first 
order equations for the so-called phase variables introduced below. 
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4.1 Computation of using phase variables 



Starting with the action (|2.1| ) we look for a black-hole solution of the form, 

ds^ = e'^^") {h-\r)dr^ + dx^^, + dt^h{r)) , (j) = (f){r) . (4.1) 

We are interested in solutions that are asymptotically AdS. In the dual field theory this 
corresponds to the presence of conformal invariance in the UV that is broken either 
explicitly by a mass deformation as in the M = 2* theory or by a marginal deformation 
as in the phenomeno logical models of |T^ . 

In the gauge 50 = one can equivalently use as the radial variable. Defining the 



following phase variables []T2| ||21|| , 

XW^'P4, i'W-il^. (4.2) 

4 V 3 A' ' AhA'' ^ ' 

the Einstein's equations can be reduced to 

This coupled first order system is sufficient to determine all of the thermodynamic 
properties (and dissipation) of the gravitational theory Once a solution to (4.3,4.4) 
is constructed, the metric functions can be determined as. 

Here 0o corresponds to the UV value at the boundary corresponding to the UV AdS 
minimum of the potential. Aq is an integration constant that essentially determines 
the energy scale of the breaking of conformal symmetry. 

The thermodynamics of the black-hole can directly be determined as follows. The 
free energy is given by 

i^(<^.) = 7^ r#/.e3^(^'')51. (4.7) 
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These backgrounds satisfy the 1st law of thermodynamics S = —dF/dT. Equation 
( [1.7| ) directly follows from integrating this equation, where 

S = J—e^^('^>^^ (4.8) 
47rG5 

is used. In the integration in ( [4.7|) one should make sure that the UV asymptotics is 
kept fixed as (ph is varied. This is explained in the case of marginal deformations in 
section 7 of 

The temperature as a function of (ph is obtained from 

T(0.) = e^(*'^) V{<P,) C . (4.9) 

iZTT 



Once we solve ( [4.3|) and ( [4.4|) above, we can calculate the free energy as a function of 
Aq and T by employing the formulae above. 

4.2 The fluctuation equation 

The fluctuation equation ( p.4|) in terms of the phase variables read, 

H';, = ci<P)H[, + d{<p)H,,, (4.10) 

where 

In passing, we note that changing the variable back to the original radial coordinate 



in ( [4.1[) produces a rather simple equation 



^„ + ^„|3^ + ^ + 2|j+i/„(^^-^|j^0. (4.13) 

where we emphasized the new terms in the bulk fluctuation eq. that arise from mixing 
of the rotationally invariant graviton excitations and the dilaton. The normalized 
frequency is defined by Uh = ourh- This equation compares with the one corresponding 
to the shear fluctuations: 

Hu + Hu i^^+l] +^12^ = 0. (4.14) 
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One crucial difference between ( 4.13 ) and ( 4.14| ) is tliat, unlike in the case of the 
shear deformation, the bulk deformation has a mass term even in the hydrodynamic 
limit Uh = 0. This implies that in general there should be a non-trivial flow from the 



horizon to the boundary in the sense of the membrane paradigm |30]. This flow is 
absent only in the case X = const which corresponds to the adiabatic limit [29|. It is 
also absent in the Chamblin-Reall solution that corresponds to constant X, see section 


In the following we apply the formalism developed here to calculate the bulk vis- 
cosity in two examples. 



4.3 Numerical results for the holographic-QCD model 

As another non-trivial example, we would like to confirm the agreement between the 
EO and the GPR formula in the improved holographic QCD model of [|l^. The model 
is based on a single scalar in the bulk theory corresponding to the operator TrF^ in 
the SU{N) gauge theory. Therefore the deformation in the UV is marginally relevant, 
hence the UV asymptotics is not of the standard asymptotically AdS type, but involve 
logarithmic corrections. In the following we present the results in the variable 

A = ev^*. (4.15) 

The scalar potential is given by, 

^(A) = {l + ^oA + V,X'/' [log (1 + V,X'/' + V,X')f'] , (4.16) 

The various parameters in ([4.16|) 

{Vo, V,, V2, V3} = {0.0413 , 14 , 5.310-9 , 170} , (4.17) 

are fixed by in order to fit the UV asymptotics of SU (N) beta-function, the observed la- 
tent heat of the confinement-deconfinement transition on the lattice and the agreement 
with the glueball spectrum in the vacuum theory 

A straightforward application of the method explained in section yields the bulk 
viscosity from the GPR formula In the figure we compare the outcome of the 
two formulae for a range of A/j that corresponds to the entire range of temperatures 
above the deconfinement transition, T > Tc on the thermodynamically favored and 
stable big black-hole branch. As seen from this figure the two formulae match to great 
accuracy. 
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Figure 2: Comparison of the bulk viscosity in the phenomenological QCD model. The 
solid (blue) curve represents the outcome of the Eling-Oz formula. The (red) dots are 
the outcome of Gubser et al's formula. 
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A 1^ in A/" = 2* plasma 

From [l^, to leading order in m^/T^ (notice the \/24 renormalization of the a to insure 
the canonical kinetic term as in (|2.1| )), 



1 

A(x) = In (53 In (2x - x^) + 5? Ai(x) , 

4 ^ ^ (A.l) 



(j){x) =V24 5i ai{x) , 

where 



ai = {2x - xY' 2^1 (|, |; 1; (1 - xf) , (A.2) 
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{2z - ^2)2 y \dy J y-1 



From ([A.l|) -( [Ar4D it is easy to deduce that near the horizon, i.e., x —> 1. 



4 /I ^, ^4 



A =ln(27rT) - + ( i - ^^1^ ) (1 - x)2 + - x)^) . 



(A.5) 



We can now compute ^ while keeping rrib fixed, see ( p.l2|) and (|2.13| ). 



B Coefficients l-Li 

In the appendix we exphcitly show the coefficients of the fluctuation equation (|375|) : 



_ (0^2(x - 1) - 12A^(2(x - 1)A' + 1))V^ (l)'\x - 1) - 3^X3 + 8{x - 1)A') 
' ~ {x - l)(f)'V '~ ^ 3{x-l)A' ' 

(B.l) 

(f)''^(x-l) -12A'(2(x-l)A' + 1) ( , , 

6(a; - 1)30'^^' V 2) 

+ 3e-2^a; V^' ] • 



C Chamblin-Reall backgrounds 

We choose an exponential potential (known also as the Chamblin-Reall geometry), p2| . 
In this case 

V = Ve^'^ ^ (C.l) 

with constant V. 

Solving the background equations we find, 

A[x) =^ - ^ , 

exp I 



6^ / + 
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where A, V are the integration constants, and without the loss of generahty we assumed 
0(0) = 0. 

To leading order in the hydrodynamic limit (|3.5|) simplifies dramatically 



= H';, + J—H[,. (C.3) 
X — 1 

We outline now the solution of the boundary value problem that we will use in the 
more complicated example of the Af = 2* gauge theory below. 

First, the general solution with the UV boundary condition (|3.6|) is given by 



H',, = 1 + K, ln(l - x) , (C.4) 

with an arbitrary constant h^v It is straightforward to rewrite ( |(J.3D in terms of 
y = 1 — X, and find the most general solution satisfying the IR boundary condition 

H'^i = hir , (C.5) 

with an arbitrary constant hi^. Matching H^^ and (the value of the function and 
the first derivatives) uniquely determines 

huv = 0, hir = l. (C.6) 

Thus, much like in 0, we conclude that for the Chamblin-Reall model 

Cn = 1 ■ (C.7) 

We may also present the results above in the language of phase variables. 
The Chamblin-Reall solution is given by (see Appendix J of 



/o 1 _ V2 



where we defined the constant. 

We note that for consistency of thermodynamics X"^ < 1, otherwise the black- hole 



solution has negative specific heat, hence corresponds to a small black-hole |21|]. One 
finds the following metric functions in the variable 0: 

m = + (c.io) 

h{(f)) = i-e'^^^-^'^K (C.ll) 
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One distinguishing fact about the above solution is that the scale function A{(j)) is 
independent of the temperature 

Before making this connection however, let us provide a simple proof — closely re- 
lated to the one given in section |C|. The fluctuation equation (|4.1CI|) simplifies drastically 
as the coefficient d{(j)) in (^4.12|) vanishes for u = 0. This means in particular that there 
is no flow from the horizon to the boundary in the sense of the membrane paradigm, 
see e.g. for the bulk- viscosity in the case of Chamblin-Reall backgrounds. 

The proof that |c]~]^|(0/i) = 1 in this case is already given in Appendix B of |29|, that 
we review here. When, the coefficient d{(f)) in ( 4.12 ) vanishes in the u = equation, 
the solution to Hn is simply given by. 



c 



dt e 



(C.12) 



<t>0 



where the function c(0) is given by ( [4.11|) and we used the boundary condition Hu (0o) = 
1. The second integration constant C is determined by the second boundary condition 
that Hii{(j)) is regular at the horizon On the other hand, the function c(0) in ( [4.11 



is positive definite because X < 0, Y > 0, < 1, and the term inside the brackets 
is given by —4/3(7, hence negative. Therefore the only way to guarantee regularity at 
the horizon is to set C = 0, hence Hu = 1 for all values of A in the limit, in particular 



-11 



is 1. 
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